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We report the first observation of tlie decay B"^ — > H^Aj witfi a significance of 8.7 a and evidence 
for tiie decay B° ~* H~A+ witfi a significance of 3.8 a. Tfie product B{B* H°A+) x B(H° 
H+tt") is measured to be (4.8lJ;9 ± 1-1 ± 1-2) x 10"^ and B{B° H^A+) x B{B~ S+tt"?!-") 
is measured to be (9.3I2 g i 1.9 ± 2.4) x 10~^. Tfie errors are statistical, systematic and the error 
of the A^ — > pK~-K'^ branching fraction, respectively. The decay B'^ — > H^A^ is the first example 
of a two-body exclusive B^ decay into two charmed baryons. The data used for this analysis 
was accumulated at the T(45) resonance, using the Belle detector at the e+e~ asymmetric-energy 
collider KEKB. The integrated luminosity of the data sample is equal to 357 fb~^, corresponding 
to 386 X 10** BB pairs. 



PACS numbers; 13.25.Hw, 14.20.Lq 

A number of iJ-meson decay modes to final states con- 
taining baryons have been observed, including h cud 
decays with either one final-state charmed meson (e.g. 

or a charmed baryon (qm. 5+ 
A~p7r"'" 0), and charmless baryonic decays that pro- 
ceed via b s OT b u transitions. Two-body baryonic 
decay modes are found to have lower branching fractions 
than multi-body modes and, in the latter, near-threshold 
enhancements are observed in the baryon-pair invariant 
mass spectra . Some theoretical models attribute these 
phenomena to baryonic form factors that are large for 
multi-body modes j^l]. 

Recently, Belle reported examples of baryonic decays 
that proceed via b —^ ccs transitions: B~ — > J/ipAp |g 
and B A+A-K 0. To date, however, nothing is 
experimentally known about two-body exclusive B de- 
cays to two charmed baryons, which would also proceed 
through b ccs transitions. An example of such a de- 
cay is 5+ 2°A+, which would proceed via the quark- 
diagram shown in Fig. ^ This two-body B decay mode, 
like B A^A^K, would produce a "wrong-sign" A+, in 
contrast to all other known B decay modes that only have 
A~'s in the final state 0. Recently the BaBar collabo- 
ration has measured the inclusive yield of (wrong-sign) 
A+'s from B decays 0. It was suggested that this type 
of B decay might be a substantial component of baryonic 
b —^ ccs transitions and could have an important influ- 
ence on the determination of the charm particle yield per 
B decay 

For exclusive two-body baryonic modes, a theoretical 
model based on QCD sum rules predicts B{B ScA+) ^ 
10~^^3- Experimental measurements of i? ^ 2cA+ test 
theoretical predictions and provide additional informa- 
tion on the dynamics of two-body baryonic B decays. 

In this Letter we report the first observation of S+ ^ 
S^A^ and evidence for B° 2^A+ decays. Charge 




FIG. 1: The quark diagram for the _B+ ~* H^Aj decay. 



conjugation is implied here and throughout the paper. 
The analysis is performed using data collected with the 
Belle detector at the KEKB asymmetric-energy e+e" col- 
lider The data sample consists of 357 fb~^ collected 
at the T(4S') resonance, which corresponds to 386 x 10^ 
BB pairs. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Chercnkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside of the coil is instrumented 
to detect mesons and to identify muons (KLM) . The 
detector is described in detail elsewhere . Two differ- 
ent inner detector configurations were used. For the first 
sample of 152 million BB pairs (Set I), a 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector were used; 
for the latter 234 million BB pairs (Set II), a 1.5 cm ra- 
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dius beampipe, a 4-layer silicon detector and a small-cell 
inner drift chamber were used We use GEANT- 

based Monte Carlo (MC) simulation to model the re- 
sponse of the detector and determine the efficiency (lif. 

We select charged pions, kaons and protons that orig- 
inate from the region dr < 1 cm, \dz\ < 4 cm, where dr 
and dz are the distances of closest approach to the inter- 
action point in the plane perpendicular to the beam axis 
(r — plane) and along the beam direction, respectively. 
Pions, kaons and protons are identified using a likelihood 
ratio method, which combines information from the TOF 
system and ACC counters with dE/dx measurements in 
the CDC 0. 

In this analysis we reconstruct the following decay 
modes: S° S~7r+ and Aif^7r+, S+ ~* S~7r+7r+, 
A^ pK^-K^ , 'Er Att^ and A pi:^ . For A — > pn^ , 
we fit the p and tt tracks to a common vertex and re- 
quire an invariant mass in a ±5MeV/c^ interval around 
the A mass. The distance between the A decay vertex 
position and interaction point (IP) in the r — cf) plane 
(dr(A)) is required to be greater than 0.05 cm and the 
angle a\, between the A momentum vector and the vec- 
tor pointing from the IP to the decay vertex, must satisfy 
cosaA > 0.995 for the case of S° AK~tt~^. We make 
no requirements on dr and |c?z| for tracks coming from 
S~ — > Att" and A — )■ p7r^ decays. For Att", we 

fit the A trajectory and the n~ track to a common ver- 
tex and require a Att^ invariant mass in a ±5 MeV/c^ 
interval around the S~ mass. We require that the dis- 
tance between the S~ decay vertex position and IP in 
the r — (j> plane to be greater than 0.01 cm. For the A's 
coming from S~ in the decay — s- 'E.~tt~^ we apply the 
requirements, dr{A) > 0.5 cm and cosaA > 0.0. For A+, 
S° and S+ we use mass windows that are ±15MeV/c^ 
around their nominal values. We use a large sample of 
inclusive A, S~, S^^° and A+ signals to verify that their 
mass peaks are well described by two Gaussians, corre- 
sponding to the core and tail of the distribution. The 
signal mass windows that are used in this analysis cor- 
respond to approximately 4ct for the core and 2a for the 
tail Gaussian. The MC studies of the inclusive A, S^, 
'E.t^^ and A+ signals show agreement with data. 

The B candidates (i.e. ScA+ combinations) are iden- 
tified by their center of mass (cm.) energy difference, 
AE = YiiEi — -Ebcam, and their beam-energy constrained 

mass, Mbc = V^^boam - i^iPi)'^^ whcrC Ebeam = ^3/2 is 

the beam energy in the cm. and pi and Ei are the three- 
momenta and energies of the B candidate's decay prod- 
ucts. We accept B candidates with Afbc > 5.2 GeV/c^ 
and \AE\ < 0.2 GeV. To suppress the continuum back- 
ground, we require the normalized Fox- Wolfram mo- 
ment [rgl R2 to be less than 0.5. We apply |cos^s| < 0.85 
for the 2° reconstruction in the AK^tt^ mode, to sup- 
press the combinatorial background. Here 6b is the polar 
angle of the _B-meson direction in the cm. 



The AE and Mbc distributions for the B+ S°A+ 
candidates are shown in Figs.|2](a) and (b), where the two 
5° modes are combined. We require Mbc > 5.272 GeV/c^ 
{\AE\ < 0.025 GeV) for the AE (Afbc) projection O. 
The hatched histograms in Figs. |21 (a) and (b) show 
the sum of normalized A+ and Sj? mass sidebands [l9| 
where no peaking structures are evident. The super- 
imposed curves are the results of a simultaneous two- 
dimensional binned maximum likelihood fit to the both 
AE versus Mbc distributions (for the two S° channels) 
with a common value of B{B^ '^c-^t) ^ ^i^c ~^ 
S+TT^) X S(A+ pK^TT^). For this fit, we constrain 
the ratio B(E° AK-tt+)/B(E° to the recent 

Belle measurement of 1.07± 0.12± 0.07 To describe 
the signal we use Gaussians with means and widths fixed 
to the values obtained from MC. The backgrounds in 
AE and Mbc are parametrized by a first-order polyno- 
mial and an ARGUS function [2^, respectively. The fit 
gives a statistical significance of 8.7cr for the signal, where 
the statistical significance is defined as a/— 21n(Lo/imax), 
where Lq and Lmax are the likelihoods with the signal 
fixed at zero and at the fitted value, respectively. The 
region AE < —0.08 GeV is excluded from the fit to 
avoid possible contributions from 5+/'^ S°A+7r"/^ 
and 5°/+ ^ SOEc^^, S°/+ ^ A+tt"/" decays, where the 
pion is undetected. The same fitting procedure applied 
separately for the two S° modes gives 12.4l3;3 (6.8 a 
significance) and 16.9l4o (5-9 o" significance) events for 
B+ S^A+ followed by 2° ^ S+tt" and B+ ^ E°^A+ 
followed by S° Aif+Tr", respectively. 

As a cross-check of the B^ ^c-^t signal, we select 
events in the B-signal region of \AE\ < 0.025 GeV and 
Afbc > 5.272 GeV/c^ for two S° modes and examine 
the A+ and mass distributions (Fig. 12 (c) and (d)). 
For the A+ (S°) distribution we require S° (A+) to be 
within ±15MeV/c^ of the nominal mass. We then fit 
each distribution with two Gaussians for the signal and a 
first-order polynomial to describe the background. The 
widths and means of the Gaussians are fixed to the values 
obtained from data as described above. The fitted signal 
yields of 32.6±7.2 events for the A+ and 29.4±6.9 events 
for the S° are in good agreement with the total signal 
yield for _B+ 2°A+, including the two S° decay modes. 

The B'^ '^c^c niode is an isospin partner of the 
5+ S^A^ mode. Therefore their branching fractions 
are expected to be of the same order of magnitude. The 
AE and Mbc distributions for the B'^ '^c-^t candi- 
dates are shown in Figs. 01 (a) and (b). The superimposed 
curves are the results of a two-dimensional binned max- 
imum likelihood fit to the AE versus Mbc distribution. 
The fit gives S.StH si gnal events. The signal significance 
is 3.8(T, taking into account the systematic uncertainty 
from the signal and background parameterization. The 
hatched histogram shows the sum of the normalized A+ 
and mass sidebands. We apply the same procedure 
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TABLE L Summary of the fit results, efficiencies, products of branching fractions and statistical significances. For the for 
two modes the product of branching fractions is B{B'^ HcAj) x S(Hc — > 3'''7r~) since for 2° AK'^tv~ we use the 
ratio B(h2 — > AA'~7r"'")/B(H° — > H~7r^) mentioned in the text. The uncertainties in the products of the branching ratios are 
statistical, systematic and the uncertainty of the A^ ^ pK^n'^ branching fraction. 



Decay Mode 


Yield 


Efiiciency(%) 


Product of B's (10^^) 


Significance 


B+^2ilA+,sil^3+7r- 


12.4+^3^ 


1.14 


5.6li:^ ± 1.1 ± 1.5 


6.8cr 






2.04 


4.0lJ:g±0.9±1.0 


5.9cr 


— > H^Aj, simultaneous fit 






4.8ii:g ± 1.1 ± 1.2 


8.7cr 


^H-AJ, =- ^H+TT-TT- 


^■^1 


0.46 


g.Slai ±1-9 ±2.4 


3.8(7 




FIG. 2: a) and b): The A_E(a)) and Mbc(b)) distributions 
for the B^ —* H^Aj candidates. The hatched histograms 
show the combined and Aj mass sidebands normalized 
to the signal region. The excess around Ai? = —0.150 GeV 
may be due to the contributions from B'^^^ — > H^A^tt"''" and 
B°/+ ^ H°E°''+, E°''+ ^ A+7r-/° decays, where the pion is 
undetected. Therefore, we exclude this region from the fit. 
c) and d): The H° (c)) and A+ (d)) mass distributions for 
the B^ —* H^A^ candidates taken from the B-signal region of 
|A£;| < 0.025 GeV and Mbc > 5.272 GeV/c^ For the H'i (A+) 
distribution we require A+ (H°) to be within ±15MeV/c2 of 
the nominal mass. The overlaid curves are the fit results (see 
the text). 



used for i?+ S^A+ to check the A+ and E~ signals 
as shown in Figs. El (c) and (d). The fit gives 9.0 ± 3.0 
events for the A+ and 8.4 ± 2.8 events for the S^. Both 
are in agreement with the i?*^ —> S~A+ signal yield. 

Table |2 summarizes the results of the fits for the 
and i?° decays, the reconstruction efficiencies including 
the B{A P7r~), statistical significance of the signals 
and extracted products of branching fractions. Here we 



FIG. 3: a) and b): The A£;(a)) and Mbc(b)) distributions for 
the — » HjAj candidates. The hatched histograms show 
the combined H~ and A^ mass sidebands normalized to the 
signal region, c) and d): The (c)) and A^ (d)) mass dis- 
tributions for the B" — > E~ candidates taken from the B- 
signal region of |A£;| < 0.025 GeV and Mbc > 5.272 GeV/c^ 
For the i-f^t) distribution we require A J (2J) to be within 
±15 MeV/c^ of the nominal mass. The overlaid curves are the 
fit results (see the text). 

use B{A+ ^ pK-n+) = (5.0 ± 1.3)% and assume 
equal fractions of charged and neutral B mesons pro- 
duced in T(4S') decays. 

The major sources of systematic error are the uncer- 
tainties in the tracking efficiency of 7% (1% per track), 
11% in charged particle identification efficiency (1% for 
pion, 2% for kaon and 3% for proton), 5% in A finding, 
6% in efficiency estimation due to MC statistics, 10% in 
the signal and background parameterization, and 13% in 
B{E° ^ AK-TT+)/B{E° E-TT+). Added in quadra- 
ture, these correspond to a total systematic error of 23% 
for B+ E°A+ and 20% for 5° ^ S-A+. 
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In summary, we report the first observation of the 
S^A^ decay mode and the first evidence for the 
^.T^^ decay mode. The products of branch- 
ing fractions B{B+ E°A+) x B{E° ^ S+tt") = 
(4.8^ J:°±l.l±1.2) X 10-5 and B{B" A+) x B{E^ 
E+TT-TT-) = {Q-^tli ± 1.9 ± 2.4) X 10-5 are measured 
with 8.7a and 3.8a significance, respectively. These re- 
sults and Belle's recent observation of i? ^ A+A-fsT 
decays are the first examples of B decays into two 
charmed baryons. The branching fraction obtained for 
'^c-^t together with the theoretic al p redictions 
for B{E° E-7r+) of - (0.9 - 2)% result in 
B{B+ S°A+) - (2.4 - 5.3) x 10"^ This can be com- 
pared with the theoretical prediction of lO"'' the 
other hand, the Belle measurement of S(J5° A+p) = 
(2.19t5]j^±0.32±0.57) x lO'^ ^ is much smaller than 
their prediction of 4 x 10"^ flU. The very large ratio of 
- 100 for BiB -> s||a+)/6(B0 -> A+p) disagrees with 
the naive expectation that the branching fractions for 
two-body baryonic B decays proceeding via b —> ccs and 
b — !■ cud transitions should be of the same order |lll | . 
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